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ABSTRACT 

We present high-contrast images of HR 4796 A taken with Subaru/HiCIAO in H-band, resolving the debris disk 
in scattered light. The application of specialized angular differential imaging methods (ADI) allows us to trace 
the inner edge of the disk with high precision, and reveals a pair of "streamers" extending radially outwards 
from the ansae. Using a simple disk model with a power-law surface brightness profile, we demonstrate that the 
observed streamers can be understood as part of the smoothly tapered outer boundary of the debris disk, which 
is most visible at the ansae. Our observations are consistent with the expected result of a narrow planetesimal 
ring being ground up in a collisional cascade, yielding dust with a wide range of grain sizes. Radiation forces 
leave large grains in the ring and push smaller grains onto elliptical, or even hyperbolic trajectories. We measure 
and characterize the disk's surface brightness profile, and confirm the previously suspected offset of the disk's 
center from the star's position along the ring's major axis. Furthermore, we present first evidence for an offset 
along the minor axis. Such offsets are commonly viewed as signposts for the presence of unseen planets within 
a disk's cavity. Our images also offer new constraints on the presence of companions down to the planetary 
mass regime (~9 M Jup at (X'5, ~3 M Jup at 1"). 

Subject headings: circumstellar matter — planetary systems — techniques: high angular resolution — stars: 
individual (HR 4796 A) 



1. INTRODUCTION 

Vega-type debris disks were first identified b y infrared (IR) 
excesses around nearby main sequence stars (Auma nn et alj 
1984). The dust content of these second-generation disk sys- 
tems is believed to be continuously re plenished via collisional 
breakup of remnant planetesimals (cf. Wyatt 2008^ 

Since the imaging of the /3 Pic system d Smith & Terrild 
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1984), nearly two dozen nearby debris disks have been 
spatially resolved. The morphological appearance of re- 
solved debris disks is predicted to be influenced by in- 
teractions between dust in the disk and near by planets 
dOzernov et al.1 120001: iKuchner & Holmanl 120031) . the local 
interstellar medium (Maness et al. 2009), rec ent stellar fly- 
by s and bi nary companio ns (Wyatt etalJl 19991) . mutual grain 
collisions (Thebault&Wu 2008), and interaction of dust 
with residual gas (iKlahr & Lira 120051) . Many resolved 
systems exhibit all morphological struc t ures predicted by 
these mechanisms dSchneider et all ! 19991: IKalas et af] l2005h 
iGolimowski et al.ll2006t iKalas et all 120081) . The observable 
morphology of resolved, optically thin debris disks is also 
wavelength dependent, as d ifferent band passes sample differ- 
ent grain size distributions dWyatJl2006l) . 

HR 4796 A is a you ng (~8-10Myr:IStaufferet alj |1995l). 
nearby (72.8 ± 1. 7 pc: Ivan Leeuwea 120071) . AOV-tvpe star 
first identified as a debris disk system from an IR excess 
observed with IRAS (Jura 1991). It has a c o-moving M-typ e 
stellar companion at a separation of 7'.'7 dJuraet al.l 119931) . 
The HR 4796 A debris disk has been spatially resolved 
at numerous optical, infrared, and sub-millimeter wave- 
lengths (e.g.. iKoerner et al]|1998l:jjavawardhana et al.l 1 1998h 
Schne ideret al|1999|lSheret et alJl200l iHinklev et al.1120091: 
Schneider et al. 2009), appearing as a narrow, highly inclined 
belt with a radius of 1'.'05. Models suggest a dual dust popu- 
lation including gravitationally confined grai ns at ~9-10 AU 
in ad dition to those imaged at ^70 AU ( Augereau et al. 
1999), p ossibly involving r adiationally evolved organic 
materials dDebes et al.1 120081) . IWahhai et all d2005l) present 
multi-wavelength observations showing that the outer disk 
cannot be modeled as a single dust component; these authors 
adopt a two-component model with two different grain sizes. 
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FIG. 1. — High-contrast images of the HR 4796 A debris disk. White plus signs mark the star's location. The saturated central area is masked in grey, (a) 
Subaru HiCIAO //-band data after reduction with simple ADI (subtraction of a median background) and median smoothing on the scale of 5 pixels (= 48 mas 
1 A/D). Orange dotted lines mark the mean of the spider diffraction pattern. The stretch is linear between ±3- 10 — 5 times the stellar peak flux, (b) The same 
after ADI reduction with the LOCI algorithm, using conservative settings (optimization region area of 10,000 PSF footprints). Note the "streamers" extending 
radially out from the ansae (yellow arrows). Linear stretch of ±1.5 ■ 10~ 5 . (c) The same with aggressive LOCI settings (300 PSF footprints) for point-source 
detection. Linear scale of ±6 • 10~ 6 . (d) HST/STIS 0.2-1.0 (im image from [ Schneider et al. 1 2009, Figure 8) reprinted for comparison. These authors point 
out a possible flux deficit in the southwestern quadrant of the disk (dashed green arrows), (e) Simple ADI image from (a) in logarithmic stretch [10~ 7 , 10~ 3 ], 
showcasing the streamers. Red shading indicates spider artifacts. 



Here we present new ground-based near-infrared imaging 
data at high spatial resolution, revealing the tapered outer 
regions of the debris disk. 

2. OBSERVATIONS AND DATA REDUCTION 

Our observations of HR 4796 A were obtained with the 
Subaru Telescope on May 24, 2011, within the SEEDS sur- 
vey (Strategic E xploration of Exoplanets and Disks with 
Subaru/HiCIAO, iTamural 120091) . The HiCIAO instrument 
dHodapp et al.l2008l) with a 20" x 20" field of view and a plate 
scale of 9. 50±0. 02 mas/pixel was used. The image rotator 
operated in pup il-tracking mode to enable angular differential 
imaging (ADI, Maro is et alJ2006l) . A sequence of 260 images 
was taken in 77-band with an exposure time of 10 s, for a total 
integration time of 43.3 min and a total field rotation of 23°. 
Weather conditions were good (seeing / /5-0 / . / 8 in I^-band) , 
and the A0188 adaptive optics system (Min owa et al.ll2.Q10h 
provided a FWHM of 6.5 pixels = 62 mas. 

The images were corrected for flatfield and field distortion 
dSuzuki et al.ll2010T) . Stellar position was estimated in each 
frame by two different methods: (1) Fitting a Moffat profile 
to the PSF halo, and (2) triangulating between symmetrical 
pairs of static speckles. Both were consistent with an empir- 
ical drift model (r$ = r Q + v i + 1/2 a i 2 ) plus measure- 
ment noise. The difference between the two sets of centroids 
showed no systematic behavior, and was consistent with in- 
coherent combination of the two measurement noise sources. 
We therefore used the drift model for image registration, for 
an estimated centering accuracy of ^0.3 pixels = 2.9 mas in 
the co-added image. 

ADI combined with the LOCI algorithm (Local l y Op- 
timized Combination of Images, Lafre niere et alj 120071) 
is currently the most successful ground-b ased imag- 
ing technique for the detection of planets dMarois et alj 
120 lOt lLagrange et"aT]l2010t ICurrie et al] 1201 lh and substel- 
lar compani ons (Thalm ann et alJ 2009; B iller et all 120101 : 
iJanson et alJ 1201 lh ! Furthermore, it has p roven useful in 
revealing high-contrast c ircumstellar disks (Thalma nn et alj 
IMotlBuenzli et al]|201Ch . 

We applied three implementations of the ADI technique to 
our data: (I) "Simple ADI", consisting of subtracting a me- 



dian background from the entire dataset before derotating and 
co-adding. This method causes the least amount of flux loss 
and is therefore useful for estimating the surface brightness 
profile of the disk. (II) "Aggressive LOCI", using frame se- 
lection criteria of 0.5 FWHM (minimum differential field ro- 
tation between images to be used for mutual subtraction, to 
limit self-subtraction of physical sources) and an optimiza- 
tion region with an area of 300 PSF footprints. Although most 
of the disk signal is lost and negative over-subtraction effects 
appear, this method achieves excellent speckle suppression 
and provides the best constraints on point sources, such as 
planets. (Ill) "Conservative LOCI", a comp romise between 
the pr evious two methods, first described in Thalma nn et alJ 
(2010). This method preserves more disk flux than aggressive 
LOCI while achieving significantly better speckle subtraction 
than simple ADI. The resulting image proves ideal for deriv- 
ing geometric parameters of the inner disk edge. Due to the 
dataset's limited field rotation budget, we do not use an in- 
cr eased frame sele c tion c riterion to conserve more disk flux, 
as iThalmann et al. (2010) do. Instead, we enlarge the opti- 
mization area to 10,000 PSF footprints, lowering the impact 
of the disk signal on the optimi zation process, as first demon- 
strated in Buenzli et al. (2010). Unless otherwise noted, we 
use the default numerical and geometric LOC I sett ings as de- 
fined in Table 1 and Fig. 1 of lLafreniere et alJ ((2007). 

3. RESULTS 
3.1. Imaging the debris disk 

Figure Q] shows the results of applying the three ADI re- 
duction methods to our data: Simple ADI in (a), conservative 
LOCI in (b), and aggressive LOCI in (c). In panel (d), we 
reprint the HST/STIS image bv lSchneider et all d2009l Figure 
8) for comparison. Starlight scattered from the disk is visible 
as a strongly projected ellipse in all images, down to stellocen- 
tric separations of 0'.'5 in simple ADI and 0'.'4 in LOCI. Dif- 
ferential imaging clearly resolves the inner edge of the disk. 

While most of the flux concentrates around this edge, each 
ansa appears radially smeared out, extending a "streamer" 
outwards along the major axis. Similar streamers have been 
observed in deep ADI or roll subtraction images of other de- 



Extended Outer Regions of the Debris Ring Around HR 4796 A 



TABLE 1 
Numerical results on HR 4796 A 



Maximum merit ellipse size 


ang. sep. 


projected (au j i ' 


Semimajor axis aP^ 


1 . Uo / ± IX UzJ 


/y.z ± 1. / 


Semiminor axis b 


u. z4o it ur. uu / 


lo.U it U.J 


Maximum merit ellipse orientation 


degrees 




Position angle £7 


26.4 zt 0.5 




Inclination i = arccos(6/a) 


76.7 zt 0.5 




Maximum merit ellipse center 


ang. sep. (mas) 


' , j /AT T\f 1 ) 

projected (AU) V ' 


x (along RA) 


— O.O ± 4.U 1 ' 


— U.4o ± U.21 v ' 


y (along Dec) 


-22.1 ±4.8< 3 > 


-1.61 ±0.28 (3) 


u (along major axis)*- 


-16.9 ±5.1 


-1.23 ±0.31 


v (along minor axis)^ 4 ' 


15.8 ± 3.6 


1.15 ±0.16 


Upper limits on planets 


5 <T contrast 


Mass (Mj up ) (5) 


At 073 = 18AUin proj. 


5.8- 10" 4 


~17 


At = 36AUin proj. 


1.2- 10" 4 


~9 


At l"=73AUin proj. 


5.7 - 10" 6 


~2.8 


At 2"= 146 AU in proj. 


8.6 - 10" 7 


~1.4 



NOTES . ( 1 ) Conversion from arcsec to projected AU assuming a distance of 
72.8 pc; the uncertainty of 1.7 p c is n ot included in the errors. (2) Based on the 
model disk described in Section l4~TI the photometric peak-to-peak radius a is 
~96% of the maximum merit semim ajor axis; thi s would yield a value of a = 
l'/044 ± 0.023 consistent with Sch neider et all <2"0"0"^) . (3) Errors in (x,y) 
are correlated; use the uncorrected errors in (u, v) instead. (4) The axes 
of the (u, v) coordinate system are aligned with the disk's major and minor 
axes, thus they are rotated counterclockwise from the RA/Dec coordinate 
system (a;, y) by H = 26. 4°. (5) Conversion from iJ-band cont rast via the 
COND evolutionary models lAllard et a l. 2001; Baraffe e t al.l2003h , assuming 
an age of 8 Myr. 



bris d isks, such as H P 61005 ("The M oth", iBuenzli et all 
120101) and Fomalhaut dKalas et alJl2010h . They are thought 
to represent part of the thin, extended, tapered outer reaches 
of the debris disk, which are most detectable at the ansae. 
Oversub traction makes them appear narrower than they are. 
In both simple ADI and LOCI, the flux at a give n separation 
is unaffected by any flux a t smaller separations (Maro is et alJ 
l2006t Lafreniere et al. 2007), thus the streamers cannot be ar- 
ti facts caused by t h e brig ht ring further inwards. 

Schneider et al.l ( 120091) note a marginal flux deficit in the 
southwestern quadrant of the debris disk. Although there are 
some indications of substructure in the same location in our 
dataset (dashed green arrows), the S/N ratio is insufficient to 
confirm the existence of a physical disturbance in the disk. 

3.2. Ring geometry from maximum regional merit fitting 

The crisp representation of the disk's inner edge in the con- 
servative LOCI image allows accurate measurement of its ge- 
ometrical properties. For that purpose, we have devised the 
maximum regi onal merit technique , which is loosely based on 
the method of IBuenzli et al.l (12010b . 

First, we apply a median filter on the spatial scale of 5 pix- 
els (psl X/D), which reduces the pixel-to-pixel noise while 
keeping larger edges and structures intact. We then convert 
the image into a signal-to-noise map (S/N; Figured) by cal- 
culating the standard deviation in concentric annuli around the 
star and dividing the pixel values in each annulus by this noise 
level. The disk itself is masked during noise calculation. The 
S/N map emphasizes the well-resolved outer parts of the disk. 

We then generate a large number of ellipses with random- 
ized major and minor semi-axes (a, b), position angles (fi), 
and center positions relative to the star (x,y). We assign a 
merit value -q to each set of parameters, defined as the mean 
value of the pixels in the S/N map that lie within an ellipse- 
shaped sampling area described by the numerical parameters. 
The sampling area for a given parameter set (a, b, Q, x, y) is 
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FIG. 2. — Derivation of geometric parameters, (a) Signal-to-noise map of 
the conservative LOCI image of HR 4796 A, calculated in concentric annuli 
around the star and excluding the bright ring from noise estimation. The 
stretch is [— 10o\ lOo - ]. (b) The same, but with red shading marking the 
sampling ellipse pixels for the best-estimate ellipse, (c) Histogram of merits 
r\ measured in empty sky, showing a quasi-Gaussian distribution, (d) Scatter 
plot of the center offset parameters (x, y) for the solution with merits r) within 
1 <j v (black), 2 a v (red), and 3 a v (orange) of the maximum merit value r] mllx , 
drawn from Monte Carlo simulations of 300,000 parameter sets. A white plus 
sign marks location of the best estimate (x, y); a black plus sign marks that 
of the star. To guide the eye, white dashed ellipses with semiaxes of 1—4 
spatial standard deviations of the 1 a v family along the major and minor axes 
are provided, as well as dotted black circles marking the 1-4 a uncertainty 
regions of the star's position. 



defined as the area between two nested ellipses that share the 
same orientation and center (O, x, y) but have semi-axes of 
(a + 5,b + e) and (a — 5, b — e), respectively. We choose 
e = 0.7px, which results in a single-pixel thickness along the 
minor axis, and S = 4px, for a 7-pixel thickness along the 
major axis. The value for 6 corresponds to the observed width 
of the high-flux areas in the ansae of the disk; this is to ensure 
that a well-fitted sampling ellipse will be delimited by strong 
gradients in the S/N image on all sides, and thus sensitive to 
displacement and scaling (Fig.|2j)). 

We find that this method yields more accurate results for 
our data than the conventional method of fitting of an el- 
lipse to flux maxima in the ring (as applied to HST data by 
ISchneider etail 12009). since our data has finer spatial sam- 
pling but a higher photometric noise level than the HST data. 
The disk profile is flat around the flux maxima and thus vul- 
nerable to noise, while the tracking of sharp gradients is ro- 
bust and profits from the fine spatial sampling. 

We find a maximum merit of r/ max = 3.46. In order to 
estimate the noise inherent to this merit figure, we measure 
r\ for 10 4 sampling ellipses placed in the empty sky region 
90° away from the position angle of the disk. We find a 
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FIG. 3 . — Modeling of the HR 4796 A debris disk, (a) Simulated scattered- 
light image of the simple power-law model disk described in the text, dis- 
played in linear stretch, (b) The same in logarithmic st retch, (c) PSF- 
subtracted HST/STIS image from Schneider et al. ( 2009) in logarithmic 
stretch. Note that a faint halo as predicted by the model is visible beyond 
the bright ring, (d) Result of conservative LOCI reduction of our Subaru 
HiCIAO data after injection of the model disk at a 90° position angle off- 
set (linear scale). The model appears qualitatively indistinguishable from the 
real disk, demonstrating that the morphology of the disk signal in our ADI 
images is well understood. 



roughly Gaussian distribution (Fig. |5J:) with a standard de- 
viation a„ = 0.077 and a mean value \i„ = —0.0006 (« 
within the errors). We therefore define the family of "well- 



a v , and the 



fitting" solutions as those satisfying r\ > rj m i 
"best estimate" (a, b, Cl, x, y) as the mean value of the well- 
fitting solution family. 

Figure[2jl plots the spread of the center offset values (x, y) 
of solutions within 1 <J n , 2a n , and 3a n of the maximum 
merit value 77 max - The contours of those solution families are 
concentric ellipses with axes aligned along the debris ring's 
axes, demonstrating approximately Gaussian error behavior. 
To further test the validity of this maximum regional merit 
tech niqu e, we produced 10 model disks (as discussed in Sec- 
tion l4.lt with random offsets from the star, injected them into 
the empty sky region in our data, reduced them with conser- 
vative LOCI, and retrieved the ellipse parameters as described 



above. Deviations between measured and real offset coordi- 
nates were consistently below 1 a, with an RMS of 0.4 a. 

The best-estimate parameters and their errors, including 
both fitting errors and systematic errors from image registra- 
tion and plate scale, are presented in Table [TJ Most notably, 
we find that the ellipse center is offset along its major axis 
(towards the southwest) by u = 16.9 ± 5.1 mas, deviating 
from the star's position at 3.3 a significance. Since the ma- 
jor axis is unaffected by foreshortening from projection along 
the line of sight, this translates directly into a physical dis- 
tance of 1.2 ± 0.3 AU , assuming an approxi mately circular 
cavity. Combined with Schneider et al. (2009)'s independent 
measurement of 19 ± 6 mas (3.2 a confidence), we can now 
confirm that HR 4796 A's debris disk is offset from its star. 
Furthermore, our data suggest an offset along the minor axis 
of 15.8 ± 3.6 mas in projection. If confirmed, this would cor- 
respond to a physical offset of ^5 AU. 

3.3. Constraints on planets 

Aggressive LOCI reduction of our image data does not 
detect any planet candidates in the system, but imposes 
upper limits on the masses of potential unseen planets. 
W e use the photometric calibration method described e.g. 
in iThalmann et al.l (1201 lb . estimating and correcting the 
separation-dep endent partial flux loss from LOCI processing 
as described in Lafreniere et al.l (120071) . and converting the H- 
ban d flux into planet mass using the COND evolutionary mod- 
els (iMEdetaDllOOll l^feetd]|^0l. Table [TJ lists the 
resulting 5 a contrasts and planet masses for four separations. 
The noise is calculated in concentric annuli, and the spider 
artifacts are masked out for this purpose. 

4. MODELING 

4.1. Simple power-law model 

In order to understand the morphology of our disk images 
and gauge their usefulness for deriving disk properties, we 
generate a model disk by calculating synthetic scattered light 
images of inclined rings, with a power-law distribution of the 
scattering cross section R. We assume the inclination of the 
disk to be given by the axis ratio, and construct a Keplerian 
disk with a small but nonzero eccentricity to match our ob- 
served semi-major axis offset (1.2AU, Table [TJ). The main 
free parameters of the disk model are the radial power-law 
slopes, 



R(r) 




-2a m 



-2a n 



-0.5 



(1) 



where r is the radial distance to the star projected onto an 
ellipse with the desired geometry, ao is the semi-major axis, 
a m > and a ou t < 0, and R(r) is related to the surface den- 
sity by cr{r) = a -R(r)-(r/ao). We use a Henyey-Greenstein 
phase function and assume it is the same throughout the disk. 
To compute images, we use the single-sca ttering radiative 
transfer code GRaTeR ( Auger eau et al J 19991) to compute syn- 
thetic images. This approach is valid for an optically thin disk. 

Given the scope of this model, we do not attempt to fit but 
merely approximate the real debris disk. Our working pa- 
rameters are a m = 35, indicating a very steep inner edge, 
a out = —10, and an asymmetry factor g = 0.10, implying a 
sl ight preference for for ward scattering (c f. q = 0.16 ± 0.0 6 
in lSchneider etld]2009t g = 0.03-0.06 in lDebes et al1 : 2008). 
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FIG. 4. — Surface brightness profiles of the HR 4976 A debris disk (black 
curves, dotted 1 a error bounds) extracted from the ADI image in Figure [5] 
in strips 1 1 pixels wide oriented along the major axis. The blue curves mark 
the corresponding profiles in the model disk in the same image. The latter 
demonstrate that the straight slope of the injected model disk is preserved in 
simple ADI reduction. The red lines are linear fits to the plots in the interval 
l'.'l— 1'/6 with respect to the disk center. The numbers indicate their slopes 
(i.e., the power-law exponents), including fitting errors. 



Figures [3^ and [3j? show the simulated scattered-light im- 
ages of the model disk at linear and logarithmic intensity 
scale, respectively. In Figure [3}:, w e show the PSF-subtrac ted 
image of the HST/STIS data from ISchneider et all (120091) in 
logarithmic scale for comparison, revealing an extended halo 
beyond the bright ring that matches the expected morphology. 
Since this dataset is reduced with PSF subtraction rather than 
ADI, it suffers from virtually no flux loss, but the background 
is dominated by residual speckle noise. 

To investigate the influence of ADI on disk flux, we inject 
the simulated scattered-light images of the model disk into 
our raw Subaru/HiCIAO dataset at a position angle offset by 
90° with respect to the real disk. Figures fjjl shows the re- 
sults of conservative LOCI applied to this modified dataset. 
The two disks appear virtually indistinguishable, demonstrat- 
ing that our simple model adequately explains the observed 
morphology. 



4.2. Surface brightness profiles 

We explore the viability of using the simple ADI image 
for the purpose of determining the .ff-band surface brightness 
profile of HR 4796 A's debris disk. We extract a strip with a 
width of 1 1 pixels centered on the disk's major axis from the 
intensity image, and collapse it along the minor axis. We re- 
peat this for the simulated disk image. The results are shown 
in a log-log plot in Figure [4] The 1 a error bounds are based 
on the radial noise profile in empty sky quadrants of the sim- 
ple ADI image before injecting the model disk. Simple ADI 
preserves the model surface brightness profile. The real disk's 
observed profile is consistent with a single power law. 

5. DISCUSSION 

We have imaged and characterized the tapered outer bound- 
ary of HR 4796 A's debris disk using ground-based angular 
differential imaging, and find that its .ff-band surface bright- 
ness profile is well described by a simple power law with a 
slope around —9.5. Our observations are consistent with the 
expected result of a narrow planetesimal ring being ground 
up in a collisional cascade, yielding dust with a wide range of 
grain sizes. Radiation forces leave large grains in the ring and 
push smaller grains onto elliptical, or even hyperbolic trajec- 
tories. Simulations of such this process predict a smoothly 
tapered outer boundary, consistent with ou r findings (e.g., 
Krivo v et alj 120061: IStrubbe & Chiang! 120061) . IWahhai et ail 
(2005) have demonstrated that the disk cannot be modeled 
with a single dust grain size, implying that dust populations 
of varying grain size must be involved. The brightness slope 
is a result of changing dust densities and grain properties with 
distance. More detailed modeling is needed to relate the ob- 
served slope to the disk's physical properties. 

Using our maximum regional merit technique, we have 
corroborated the evidence for an offset between the debris 
disk's inner edge and the star. The dynamical influence 
of unseen planets orbiting within the disk cav ity is most 
commonly invoked to expla i n such offsets (e.g., iKalas et alj 
120051: iThalmann etafl 120091: iBuenzli et all 120101 and refer- 
ences therein), and may be the cause of the proposed but un- 
confirmed intra-ring gap. Our confirmation of the ring offset 
also speaks against a dyn amically cold source p lanetesimal 
population as proposed by Thebau lt & Wul (120081) as an alter- 
native to the planet hypothesis for explaining the morphology 
of the ring. The planets may well be detectable with the up- 
coming next-generation high-contrast imaging facilities. 
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